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Abstract. Using first-principles density functional theory in the implementation of 
the exact muffin-tin orbitals method and the coherent potential approximation, we 
studied the surface energy and the surface stress of the thermodynamically most stable 

' ^ ■ surface facet (100) of the homogeneous disordered body-centred cubic iron-chromium 

system in the concentration interval up to 20 at.% Cr. For the low-index surface facets 
C_> ' of Fe and Cr, the surface energy of Cr is slightly larger than the one of Fe, while the 

surface stress of Cr is considerably smaller than the one of Fe. We find that Cr addition 

to Fe generally increases the surface energy of the Fe-Cr alloy, however, an increase of 

^ ' the bulk amount of Cr also increases the surface stress. As a result of this unexpected 

trend, the (100) surface of Fe-Cr becomes more stable against reconstruction with 

rvQ ' increasing Cr concentration. We show that the observed trends are of magnetic origin. 

In addition to the homogeneous alloy case, we also investigated the impact of surface 

■^ ' segregation on both surface parameters. 
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1. Introduction 

In recent years, iron rich Fe-Cr alloys, as the basis of ferritic (body-centred cubic (bcc)- 
based) stainless steel, have attracted much scientific attention for their potential use in 
the next generation fission and prospective fusion reactors. [1] Employed as first wall 
and blanket material, this steel must withstand neutron-induced radiation damage, [2] 
e.g., swelling and void formation. Moderate Cr addition to bcc Fe in the order of 
10 % most beneficially improves its swelling and irradiation creep behaviour. [Sj H] 
Chromium substitution is further known to influence other mechanical properties such 
as the ductile-brittle transition temperature and radiation-induced hardening. [5] The 
experimental corrosion resistance of ferritic steel in an oxidising environment improves 
drastically if the bulk alloy concentration exceeds 9-13 weight percent Cr and this self- 
healing protection attributed to the formation of a chromium oxide scale. [6l [7] There 
is, however, a limit to the amount of Cr that may be added to steels as the beneficial 
low corrosion rate is shadowed by the enhanced precipitation of intermetallic phases 
which often degrade the mechanical properties of ferritic steel. Moreover, it has also 
been recognised that alloying with additional elements further increases the pitting and 
crevice corrosion resistance in certain aggressive environments, e.g., molybdenum in 
chloride environments. [8] 

The surface of the bcc Fe-Cr system has often been used as the prototype reference 
system to study the behaviour of surfaces of ferritic stainless steels. There is a great 
deal of phenomenological modeling on the passivity of stainless steels, [9l [TOl [11] while 
first principles investigations of surfaces focused on the atomic level behaviour involving 
surface segregation and the surface magnetic structure for the technologically relevant 
Fe-rich Fe-Cr alloys. [T2l[l3l[Tl[T5l[T6l[I71[l8l[T9l|20] Understanding the aforementioned 
threshold behaviour of the corrosion resistance and the particular role of Cr in the 
segregation process has been a primary target for modeling. Experimental evidence 
indicates Cr enrichment at the surface at high temperatures and bulk Cr concentrations 
larger than 13at.%. [211 l22l [231 121] Recent ab initio calculations indeed predicted a 
sharp transition from Cr-free to Cr enriched surfaces at around 9at.% Cr in the bulk 
alloy, [ini HH [12] attributed to complex magnetic interactions between ferromagnetically 
interacting Fe and anti-ferromagnetically interacting Cr species, which are likely to 
govern many essential characteristics of the Fe-Cr phase diagram below the Curie 
temperature. [25] 

Not much is known about the ferritic Fe-Cr system concerning two essential 
macroscopic parameters that describe the thermodynamic properties of its crystalline 
surface: surface energy and surface stress. That is surprising since the significance of 
stress and strain effects on surface physics has been widely discussed. [26], |2I1 [23 [22] 
The equilibrium shape of mesoscopic crystals is the one that minimises its surface free 
energy. [30] The surface energy is further of eminent relevance in connection to faceting. 
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roughening, crystal growth phenomena at surfaces and has been discussed in relation 
to segregation. There has been an increasing experimental and theoretical activity 
to understand the importance of stress on many physical properties associated with 
surface relaxation and reconstruction, [311 [32] segregation, [33] surface adsorption, [31] 
and its role in bottom-up self-organisation and surface melting. [26] On the other hand, a 
theoretical study of surface parameters of crystalline surfaces of the individual elements, 
Fe and Cr, was subject to a number of publications. [351 ESI |371 |3H1 133 According to 
the general expectation, surface stress of clean surfaces is tensile due to the increased 
electron density within the surface layer. However, it was demonstrated recently [361137] 
that magnetism can overwrite this picture leading to exceptional, compressive surface 
stresses as predicted for the thermodynamically stable surfaces of bcc Cr and cubic Mn, 
though not in the case of bcc Fe. It is hence worthwhile to investigate the surface of 
the bcc Fe-Cr system to gain information on its essential surface parameters keeping 
in mind the presence of complex magnetic interactions in this alloy system that may 
alter our expectation on their behaviours. In this context, it is important to mention 
the particular role of ab initio calculations in the determination of surface parameters 
since experimental methods to determine their absolute value often lack reliability and 
accuracy. [271[28lll0lllllll2lll3lllllll5] 

This work deals with an ab initio determination of surface parameters for the 
thermodynamically stable surface of bcc Fe-Cr ((100) facet) in the concentration range 
of 0-20 at. % Cr. The paper is organised a follows: In Section |2] we briefly overview 
the theory of surface energy and surface stress. Numerical details of our computation 
are presented in Section [31 We discuss our results in Section H] for two different surface 
chemistries: a perfectly truncated bulk system without spacial concentration dependence 
and a system involving surface segregation. The reason for the latter is to account for 
the observed transition from Cr-free to Cr enriched surfaces around the aforementioned 
threshold bulk Cr concentration. 

2. Surface parameters 

Surface energy and surface stress are two fundamental quantities to characterise the 
macroscopic properties of surfaces. Qualitatively, the scalar surface free energy, 7, was 
introduced by Gibbs as the reversible work per unit area to create a surface. [271 EH] The 
tensorial surface stress, Tjj, i,j = {x,y}, is the reversible work per unit area to stretch 
a surface elastically in the surface plane which is here assumed to lie in the x-y-plane. 
If 7 < for a particular surface of a solid, then this surface is unstable and the crystal 
fragments spontaneously. Hence 7 is positive for stable bulk systems. The components 
of T may be both positive (tensile surface stress) or negative (compressive surface stress). 
Tensile surface stress on a surface favours smaller in-plane lattice constants than the 
bulk value while a surface with compressive surface stress favours a larger one. 



Surface parameters of ferritic iron-rich Fe-Cr alloy 

Both fundamental quantities are related by Shuttleworth's equation 

Q'-y 

Tij = l5ij + -r7-, (1) 

where eij denotes the components of the stress tensor specifying an elastic in-plane 
deformation of the surface {5ij Kronecker delta). The absolute value of the residual 
surface stress, d^jde^j^ may exceed the specific surface energy. [271 l2H] Its is obvious 
that negative surface stress in solid materials must hence result from a negative residual 
surface stress (whose absolute value surpasses 7 for the diagonal components). 

In ab initio total energy calculations, the surface energy is usually computed as, 

7 = ^ , (2) 

where i?surf and -Ebuik specify the energy of two semi-infinite bulk systems and the 
infinite bulk system, respectively, normalised to the unit area A. Surface energies are 
conveniently extracted from slab calculations and different procedures were proposed to 
yield convergent numbers with the slab size [ITJ HH HO] . In the present work, we model 
the Fe-Cr system by considering two distinct subsystems: one that includes the surface 
(surface subsystem) and one without (bulk subsystem). Due to the periodic boundary 
conditions parallel to the surface, the size of the slab has to be converged with respect to 
the thickness of the slab only. Here, we follow essentially [37] and derive both -Esurf and 
-^buik from slabs with the same thickness characterised by the total number of layers, 
ra, taken as a multiple of the bulk equilibrium lattice parameter oriented normal to the 
surface plane. In case of the surface subsystem, the slab consists of an atomic part with 
thickness n^ and vacuum which is needed to decouple the two surfaces of the slab from 
another (across the vacuum). This surface-surface distance is denoted by n^. Since the 
slab representing the surface subsystem contains two equal surfaces, a factor of one half 
is added to ([2]) to yield the surface energy of one surface, viz. 

^ = 2A ' (') 

where n = n^ + n^, and both -E"m,f and El^^^^ refer to the total energy of the n-layer 
slab. -Ebuik is scaled by a factor of n^/n to the correct number of atoms. Only geometry 
relaxation in the direction perpendicular to the surface is allowed and may be included 
in -E^rf as appropriate. 

One out of many equivalent definitions of surface stress is to define it as an excess 
stress of the bulk stress tensor, a}j, at the material- vacuum interface, [261 EH [28] viz. 

r-O 



Tij 



j(^J^dAjjz[a,,{z)-a^y (4) 

where the surface excess quantity is defined per unit area A. The components of the 
surface stress tensor, crjj(z), depend on the coordinate z oriented perpendicular to the 
surface and indicate the change of the bulk counterpart in the semi-infinite space ranging 
from z : [—00 : 0] due to the existence of a surface at 2; = 0. The in-plane integration 
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dA = dxdy extends over the surface area A. aij{z) and o"|^are assumed hereafter to have 
no lateral spacial dependence. 

The change of the elastic energy out of equilibrium due to an infinitesimal 
deformation 5eij of the system that possesses the surface is [26l [50] 



^El^ri = A I dz ^a,j{z)5eij, (5) 



u 



which may be rewritten with the help of (HI), 

The second term on the right hand side is however just the elastic energy of the 
bulk system, SEJ^^^y^. Differentiating the previous equation with respect to Cjj at the 
unstrained state yields 

1 ^{Elnri - ^bulk) 



'''' ~ A de, 



(7) 

For a high-symmetry surface facet such as the bcc (100) facet, t^x = Tyy, r^y = Ty^, and 
if we further assume an isotropic distortion, e^x = ^yy = ^ (zero otherwise), we arrive at 

^ _ ^ 1 ^(^s'urf - ^bulk) 

r = Tii — 



2A de 



(8) 

=0 

A factor of 1/2 appears due to the apphed isotropic strain. Equation ([7]) (or its more 
specific form, (jS])) is conveniently [501 ESI EZl EB IS2] used to compute the surface stress 
from the elastic energies of a surface subsystem and a bulk reference system employing 
slabs. As for the computation of surface energies, -E^m-f and E^^^^ are in fact computed 
for the n-layer slab, and the bulk energy is scaled by a factor n^/n to the number of 
atoms in the surface subsystem. Also, another factor of 1/2 is introduced in ([8]) due 
to the presence of two surfaces in the surface reference slab. Numerically, the elastic 
energies in ([8]) are fitted to second order polynomial functions with fit coefficients c and 
rf as a function of strain (e), 

E,^^f(e) - E^urf (0) = Csurfe + rfsurfe' (9a) 

71 

-^ (El'-^ie) - EUm = Cbuike + 4uike^. (96) 

The reference energies, -^^rf /bulk ('-')' ^^^ ^^^ energies of the unstrained surface and 
bulk states, which may serve to compute 7 following ([3]). Assuming the previous fit 
functions, we may express r by 

Csurf Cbulk /-in\ 

"= AA • ^'°) 

Although in theory Cbuik = for bulk in equilibrium, it may be finite in calculations due 
to numerical errors. [37] 

We conclude this section by mentioning the excess surface stress, [27l [28l [53] which 
may be derived from ([T|). For diagonal components, 

|-Hn.-,|. (11) 
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A larger value of the excess surface stress, {th — 7I, indicates a higher tendency of a 
surface towards reconstruction. [53l [5i 



3. Electronic structure calculations 

Total energy calculations within the framework of density functional theory (DFT) 
were done by means of the exact muffin tin orbitals (EMTO) method, [55l [561 157] 
which is a screened Korringa-Kohn-Rostoker type of method [581 ESI EO] and solves the 
Kohn-Sham equations in a Green's function formalism. This enables us to compute the 
electronic structure of substitutionally disordered alloys using the coherent-potential 
approximation (CPA) [611 ESI ES]. The one-electron potential is represented by 
optimised overlapping muffin-tin potential spheres describing more accurately the exact 
crystal-potential compared to conventional muffin-tins or non-overlapping spherical 
symmetry potentials. [551 El] Combined with the full-charge density technique [651 EQ 
[671 168] for total energy calculations, [551 ISH] the EMTO method has proven to yield 
reliable total energies and electronic structure in practise including the Fe-Cr system 
under consideration. [69l[ll[l3lI70l[7IlI72l[73lI711ll5] The CPA, being a single-cite 
approximation to the impurity problem, is a standard technique for electronic structure 
calculations in totally random alloys, suited for the case of alloy components having 
similar sizes. Due to its single-site nature, screening corrections to the potential and the 
total energy must be taken into account which was done within the screened impurity 
model of Korzhavyi et al [76] 177] 178] The CPA cannot directly treat atomic short range 
order (ASRO), which was observed for the Fe-Cr system by means of neutron diffusive 
scattering. [79l |80] Accordingly, there is a tendency to form an ordered compound for 
Cr concentrations smaller than llat.% while larger Cr contents incline short-range 
clustering. Inclusion of ASRO is, however, beyond the scope of this work. 

Total energies in EMTO were obtained using the Perdew-Burke-Ernzerhof (PBE) 
parametrisation [8T1 [82] of the exchange-correlation energy functional, while self- 
consistent charge densities were computed in the local-density approximation in the 
parametrisation of p3]. Switching off gradient corrections to the exchange-correlation 
potential is justified since functionals in the generalised-gradient approximation are 
known to overestimate the magnetic moment of iron. This perturbative approach gives 
accurate total energies which was also tested in the case of Fe-Cr. [75] [8^ The EMTO 
partial waves were expanded into s, p, d, and / orbitals, and the core states were 
recalculated at each iteration step. Integration over the Brillouin zone was done on a 
fc-point grid of 15 x 15 x 1 points in case of the surface subsystem. A single fc-point along 
the short reciprocal lattice vector (corresponding to the direct lattice vector parallel to 
z) is sufficient since bands are dispersionless in this direction (ensured by the converged 
thickness of vacuum). For the bulk reference system, the number of fc-points in this 
direction was converged to the value of two. Increasing the /c-point grid of the surface 
subsystem (bulk subsystem) to 20 x 20 x 1 (20 x 20 x 3) showed that the total energies are 
converged at a level of 1 meV/atom. The Green's function was evaluated for 16 complex 
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points distributed exponentially on a semicircle including states below the Fermi level. 

For the particular case of Fe we also carried out DFT calculations with the 
full-potential local-orbital scheme, FPLO-9. 01-35 |H5] and PBE. The convergence 
of numerical parameters and the basis was carefully checked. Linear-tetrahedron 
integrations with Blochl corrections on a 12 x 12 x 1 (12 x 12 x 3) mesh in the full 
Brillouin zone for the surface subsystem (bulk subsystem) were sufficient to converge 
the total energy at a level of 1 meV/atom compared to a 18 x 18 x 1 (18 x 18 x 5) mesh. 
The valence basis of Fe comprised 3spd, Aspd, and 5s states. 

Ferritic Fe-Cr crystallises in the bcc crystal structure, and the (100) surface was 
reported to be the most stable facet. [I3] In case of EMTO, we modelled this surface by a 
slab with a converged thickness of 13 atomic layers (approximately 37 A) and separated 
by vacuum with a thickness oi n^ = 7 (approximately 20 A). A somewhat larger surface 
subsystem with n^ = 19 and n^ = 7 was required to yield converged surface parameters 
in case of FPLO. The symmetry of the slabs includes a mirror plane parallel to the 
(identical) surfaces. All atomic positions were held fixed to the ideal bcc lattice sites. 
Surface geometries may differ from ideally truncated bulk crystals since relaxation and 
reconstruction may occur. Surface reconstruction on the close-packed surfaces of the 
transition metals is rather uncommon, [86] however relaxation of the surface layer and of 
subsurface layers is frequently observed. Results of [87] and references therein indicate 
that layer relaxations of the most stable (110) and the second most stable (100) surface 
facets of Fe are minor and change the corresponding surface energies in the order of 
1% (cf. comparison in table [2]). Punkkinen et al recently compared the surface energy 
and the surface stress of the most stable surfaces facets of Fe and Cr for non-relaxed 
geometry with values of fully relaxed surface geometries. [36] Due to the enhanced 
surface magnetism in both systems, relaxations are only minor and they were found not 
to markedly alter the surface energy and the surface stress. We expect for the same 
reasons that the surface geometry of the bcc (100) surface of Fe-Cr remains close to 
the truncated bulk one. To further support this point, we exemplary computed the 
surface layer relaxation of the bcc (100) surface of Fe with EMTO. We found a relative 
change in length of the interlayer distance due to relaxation of —2.4% (negative sign 
indicates inward relaxation) and an accompanied reduction of the surface energy by 
1.1 %. The available experimental data as obtained from low-energy electron diffraction 
on the top layer relaxation of the bcc (100) surface of Fe is ambiguous with values of 
-1-0.5 % and —1.4 ± 3%, both from [HE], and —5 ± 2% from [89]. One has, however, to 
keep in mind that our computation was restricted to the relaxation of the top-layer only 
while subsurface interlayer distances remained rigid. The total energies as a function of 
strain, cf. ( l96l) . were computed in a strain interval of |e| = 0.02. 

We conclude this section by establishing the precision of the calculated surface 
parameters in EMTO with respect to the selected size of the surface subsystem 
{n^ = 13) and the strain interval (|e| = 0.02) using the example of Fe. The previous 
set of numerical parameters define the reference system, which we consider to yield 
converged surface related quantities in this work. The surface energy and the surface 
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Table 1. Influence of the choice of the size of the surface subsystem, characterised by 
n,„, and the strain interval for elastic energy fits, |e|. The surface parameters of the 
reference system, given by rim = 13 and |e| = 0.02, are highlighted in boldface. The 
absolute difference A as well as the relative change in % in parenthesis are specified, 
and stated at the end of each row and at the end of each column for a decrease of 
|e| from 0.02 to 0.01 and an increase of rim from 13 to 15, respectively. 7, r, and the 
absolute difference are in units of J • m~^. In all cases n^ = 7. 



rim 


surface energy 7 


surface stress r 


A 




|e| =0.02 


|e| =0.01 




13 


2.615 


0.57 


0.50 


-0.07 (-12) 


15 


2.626 


0.51 


0.56 


0.05 (10) 


A 


0.011 (0.4) 


-0.06 (-11) 


0.06 (12) 


- 



stress of Fe were computed for a larger surface subsystem {rim = 15) on the one hand, and 
the surface stress was fitted to total energies from a narrower strain interval (|e| = 0.01) 
on the other hand. Table [1] lists both the absolute values of the surface parameters 
and the change (A) with respect to the reference set of surface parameters. Apart from 
pure Fe, we assessed the precision for several alloy concentrations of Fe-Cr in the same 
way, and noted that the tabulated A's in table [H represent characteristic values for all 
concentrations tested. A may be used to define the precision of our calculation with 
respect to the choices of Um and of |e|, which is hence of the order of 0.01-0.02 J ■ m~^ 
for surface energies, and approximately 0.04-0.08 J- m~^ for surface stresses. The energy 
scale for surface stress calculations is roughly one order of magnitude smaller than the 
one for surface energy calculations which explains the difference in the corresponding 
A's. 

Due to the neglect of ASRO in this work, we model a chemically homogeneous bulk 
alloy, i.e., there is no spatial probability (composition) dependence of the distribution 
of the alloys components. First we consider an ideally truncated bulk system, i.e., a 
system with surface, for which no atomic redistribution (segregation) occurs. Hence the 
surface composition is identical to the one in the bulk for all compositions. Second, we 
allow for a change of the surface chemistry. 

4. Results and discussion 

4-1- Surface parameters of Fe 

For the theoretical equilibrium lattice parameter of ferromagnetic (FM) bcc iron, we 
obtained 2.837 A and for the spin moment a value of 2.21 /^b- This is to be compared 
to experimental values, 2.867 A |90] and 2.21 /iB for the total magnetic moment. [9T] 

Table [2] lists the surface energy and the surface stress of the bcc (100) surface of 
FM iron as calculated in this work and compared to available data from the literature. 
The multitude of comparable ab initio data from full-potential (FP) and projector- 
augmented-wave (PAW) methods allows to draw conclusions on typical scatter of surface 
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Table 2. Surface parameters of the bee (100) surface facet of FM iron. All methods 
employed gradient corrected density functionals. 



method 


surface energy 7 


surface stress r 


references 




(J.m-2) 


(J • m- 


-') 




EMTO, PBE 


2.62, 2.59" 


0.57 




this work 


FPLO, PBE 


3.09, 3.07" 


1.15 




this work 


PAW, PBE 


2.55^ 2.50" 


1.39" 




m 


FP-LAPW, PBE 


2.6^ 


- 




m 


PAW, GGA[92l Hg 


2.48, 2.47" 


- 




m 


PAW, GGA[92l |93] 


2.32, 2.29" 


- 




m 


experiment 


2.41^= 


- 




M 


° surface layer relaxation included 








^'estimated from figure 








■^estimated at T = K 









energy and surface stress calculations, as well as a critical evaluation of the present 
results. Concerning the surface energy of Fe, the only outlier seems to be the value 
obtained with FPLO, since the remaining surface energies scatter in the range of 
approximately 2.3-2.6 J ■ m~^. The particular choice of the gradient corrected density 
functional may have an effect on the surface energy as all PBE values are larger than the 
values obtained with the parametrisation of Perdew et al |92l |93]. Our EMTO value 
of 7 = 2.62 J ■ m"^ (without surface layer relaxation) is in close agreement to the PAW 
and the FP linear augmented plane wave (FP-LAPW) results of Punkkinen et al The 
too high surface energy from FPLO may be related to a too strongly contracted wave 
function at the bulk- vacuum interface, cf., e.g., the analysis in [Mj. Concerning the 
surface stress of bcc Fe, the FP and PAW values scatters in the range of approximately 
1.1-1.4 J ■ m~^. The EMTO value, r = 0.57 J ■ m~^, is comparatively small and may 
thus indicate a systematic underestimation of the surface stress in Fe and the Fe-based 
system. This underestimation may be ascribed to the muffin-tin approximation to the 
one-electron potential. 

The experimental value of the surface energy of Fe from [95], 7 = 2.41J ■ m~^, 
is estimated from surface stress measurements of the liquid-vacuum interface at the 
melting temperature and surface stress measurements of the liquid-solid interface, and 
extrapolated to T = K. 

Because of the reduced coordination number, the magnetic moment of Fe at the 
surface is enhanced compared to the bulk. We obtained a surface magnetic moment of 
2.97 /iB in very close agreement with a recently reported value, 2.96 /iB- [12] Calculations 
of Punkkinen et al [361 137] suggested an almost linear relationship between the magnetic 
moment enhancement, Am?, and the magnetic surface stress, Tmag, on the basis of their 
computed surface stresses for the most stable surfaces of magnetically ordered Cr, Mn, 
Fe, Co, and Ni, i.e., 

r^ag oc Am^ = m^^^ - m^^i^, (12) 
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where rrisurf and rribuik are the magnetic moments at the surface and in the bulk, 
respectively. The magnetic contribution to r, r^ag, is defined as the difference 
between the nonmagnetic and the magnetic values of r, that is evaluated without spin- 
polarisation (Tnsp) and with spin-polarisation (rgp) for identical surface geometry, viz. 
Tmag = Tnsp " ^sp- The geometry of the spin-polarised system is the reference state if not 
stated otherwise. The above proportionality was verified for elements with FM order 
(Fe, Co, and Ni) and anti- ferromagnetic (AFM) order (Cr and Mn). The present EMTO 
values for FM Fe are Am"^ = 3.91 /ig and Tmag = 3.14 J -m"^, respectively, which fit very 
well to the correlation established by Punkkinen et al (see figure 3 from [36j). 

4-2. Chemically homogeneous Fe-Cr alloy 

First we consider the case of chemically homogeneous surface alloys, i.e., it is assumed 
that the chemical composition at the surface is identical to the bulk composition. 

4.2.1. Lattice constants and surface parameters The theoretical equilibrium lattice 
parameters of ferritic Fe-rich Fe-Cr alloys (0-20 at. % Cr) were previously calculated 
with EMTO-CPA and discussed in detail in [721 [731 UH UHl E]. Since our computed 
lattice parameters practically reproduce these earlier results, we refer the reader to those 
references. It is, however, important to point out the non-linear behaviour of the lattice 
parameter of the Fe-Cr system. 

In the atomic concentration range of 0-20 at. % Cr, we identify a clear trend of the 
concentration dependence of all surface parameters in Fe-Cr alloys, see figure [H The 
surface energy increases monotonically by 0.33 J-m~^ for an increase of the concentration 
of Cr from to 20at.%. This trend is not entirely unexpected since the surface energy 
of the non-relaxed bcc (100) surface facet of Cr was found to be larger by 0.5-0.8 J • m~^ 
in theory than the one of Fe. [37] That is, the surface energy of the disordered alloy 
with low Cr content (0-20 at. % Cr) follows a monotonic trend (rule of mixing) given by 
the boundary values of pure Fe and pure Cr (the exact PAW values from [37] for the 
relaxed surface are 2.50 J ■ m~^ for Fe and 3.06 J ■ m~^ for Cr). 

In previous theoretical considerations for the most stable surfaces [361 ST], the 
surface stress of Cr was reported to be 1.9 J ■ m~^ smaller than the corresponding value 
of Fe for relaxed surface geometries and likewise was the difference for non-relaxed 
geometries. As further stated, this difference amounts to 1.7 J ■ m~^ if specifically the 
relaxed bcc (100) surface facet is considered (the exact PAW values from [37] for the 
relaxed surface are 1.39 J-m~^ for Fe and —0.32 J-m~^ for Cr). Since relaxation seems to 
have a similar effect on r for both elements for the most stable surface, it is reasonable 
to assume that the effect of relaxation on r for the (100) surfaces of Fe and of Cr are 
also similar. Hence, the surface stress of the non-relaxed bcc (100) surface of Cr is 
presumably still considerable smaller than the one of Fe (roughly by 1-2 J ■ m^^). On 
this ground, we expect an overall decrease of the surface stress of Fe-Cr with increasing 
Cr content. Our findings in the dilute Cr concentration range are however contrary 
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Figure 1. Surface parameters (in units of J/ni^) of Fe-ricli Fe-Cr alloys (0-20 at.% 
Cr): (top panel) surface energy 7, (middle panel) surface stress r, (bottom panel) 
excess surface stress |t — 7I. Lines are a guide to the eye. 
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to this expectation (figured]). In the range up to lOat.% Cr in the iron matrix, the 
surface stress strongly increases to a maximum value of 1.78 J-m~^ being approximately 
1.25 J ■ m~^ larger than the corresponding value of Fe. r levels off for concentrations 
higher than lOat.% Cr. 

The third surface-characteristic quantity depicted in figure [1] is the excess surface 
stress, |r — 7|. It evidently decreases strongly in the concentration range up to lOat.% 
Cr relative to the value of pure Fe, which is mainly due to the accompanied increase 
of T. For Cr concentrations above lOat.%, \t — 7I remains almost unchanged. The 
surface reconstruction is predicted to occur when the excess surface stress becomes larger 
than the characteristic surface strain energy associated with the reconstruction. [531 EU 
The latter may be expressed in terms of the shear modulus and the Burgers vector. 
Now, taking into account that the elastic moduli of Fe-Cr alloys show a rather weak 
composition dependence for the present concentration interval, ^74j one may conclude, 
that the (100) surface of Fe-rich Fe-Cr alloys is considerable more stable against 
reconstruction than the surface of pure Fe. 

4-2.2. Magnetism and magnetic surface stress The magnetic structure of ferritic Fe- 
Cr is governed by interactions between Fe atoms, that prefer to align their magnetic 
moments in parallel, and Cr atoms, that favour an anti-parallel alignment. Its energetics 
is rather well described within collinear magnetism of fixed Ising spins. [251 EH [75] In the 
iron-rich ferritic Fe-Cr alloys, the moment at Cr sites are coupled anti-parallel to the ones 
of Fe necessarily implying that they are aligned in unfavourable parallel orientation with 
respect to other Cr atoms. In the dilute limit, however, the average Cr-Cr distance is 
large and their mutual interaction energy small. [HI [15] The EMTO-CPA spin moment 
at Cr is — 1.62/iB on the impurity level. This calculation was done with 0.05 at. % Cr. A 
negative sign of the magnetic moment indicates an antiparallel alignment with respect 
to the moment of Fe which was defined to possess a positive sign. Klaver et al obtained 
a slightly larger spin moment of Cr in the dilute limit, — 1.8 /ib, calculated with the PAW 
method and PBE for a super cell with an effective Cr concentration of 0.19 at. %. [15] 
Thus, a Cr atom in a surrounding Fe matrix at very low Cr concentrations is much 
stronger polarised than in AFM ordered pure Cr (the experimental spin moment of the 
long wave spin-density ground state of Cr is approximately 0.59 /xb)- 

Increasing the Cr concentration beyond dilute levels results in a gradual loss of the 
modulus of the magnetic moment at Cr sites while the Fe spin moment hardly changes, 
see figure [21 These findings are in line with previously published assessments. [73l[T5l [70] 
This concentration dependent effect on the Cr moments seems to be well understood 
on the basis of an increased number of unfavourable Cr-Cr interactions with increasing 
number of Cr atoms in the Fe matrix. The total net magnetic moment of the alloys 
decreases in the same concentration interval. 

Both the magnetic moments of Fe atoms and of Cr atoms located at the surface 
are enhanced with respect to their bulk values. The EMTO-CPA magnetic moment 
of Cr located at the surface in the dilute limit amounts to — 3.11/xb being thus even 
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Figure 2. Magnetism in Fe-Cr alloys for the surface and bulk reference systems. 
Top panel: site resolved spin magnetic moments in the bulk and at the surface; 
bottom panel: concentration averaged spin magnetic moments in the bulk and at the 
surface, and the difference between surface and bulk magnetic moments. EMTO-CPA 
calculations on the impurity level were done for 0.05 at.% Cr. Lines are a guide to the 
eye. 
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larger in absolute value than the corresponding value of Fe (2.96 /ie)- We realise from 
figure [2] that Fe moments in the surface layer and in the bulk change little as a function 
of concentration in contrast to the magnetism at the Cr sites. The modulus of the spin 
moment of a Cr atom localised at the surface undergoes a slight decrease in the range of 
increasing Cr content from to 20 at. %, which is in fact similar to the reduction of the Fe 
surface moment. The different behaviours of the Cr-Cr interaction is due to the different 
average distance between two Cr atoms at the surface and in the bulk (for the same 
concentration it is larger at the surface) and the interaction energy which scales with 
the number of atoms in nearest neighbour shells (which is larger in the bulk). [14J As 
discussed above, the absolute value of the bulk Cr moment drops considerably resulting 
in a drastically higher moment enhancement of Cr at the surface. This strong moment 
enhancement is in fact a propensity of atoms in the surface layer only; Fe and Cr 
magnetic moments in subsurface layers of the bcc (100) surface possess almost bulk 
values. [T1[I2] 

As a consequence of the distinct moment behaviours of the individual alloys 
components, the total net surface moment of the alloy diminishes more pronounced than 
the total net bulk moment (figure [2]). Their difference drops to zero at approximately 
20at.% Cr, i.e., the net surface moments equals the net bulk moment. We suggest that 
the changes in the magnetic structure determine the trends of 7 and r as we argue 
below. 

To understand the contribution of magnetism to the noticed behaviour of surface 
parameters we return to figure [1] where we included data of non-spin-polarised 
calculations as well. These were done for exactly the same geometry as the spin-polarised 
calculations. The magnetic contribution to the surface energy, 7mag, is likewise defined 
to Tmag as the difference between the nonmagnetic and the magnetic values of 7. In 
agreement with previous investigations for 3d transition metals, [361 ES] we find that 
magnetism generally reduces surface energies and surface stresses in the case of Fe-Cr. 
The magnitude of the magnetic contribution to the stress is clearly larger than its effect 
on the surface energy. Furthermore, Xmag reduces strongly by approximately 1.3 J ■ m~^ 
in the concentration range from to lOat.% Cr while it is approximately constant 
for higher concentrations. 7mag exhibits the same behaviour, however less pronounced. 
These behaviours confirm that changes in the magnetic structure drive the observed 
trends of 7 and r. 

Both Fe and Cr are polarised in Fe-Cr, hence both species contribute to r^ag- To 
better understand their individual contributions, we explore the concentration weighted 
surface moment enhancement, that is, we evaluate ( fT2l) separately for Fe and Cr and 
weight the results by (1 — x) and x, respectively. The resulting data in figure [3] signals 
a correlation between (1 — x)Am'^{Fe) and r^ag in the concentration range x < 20 at.% 
Cr. The effect of Cr on Tmag is strongly diminished since the monotonically increasing 
weighted surface moment enhancement of Cr, a;Am^(Cr), is not reflected in the trend 
of Tmag- Chromium may, however, be associated with the levelling-off of the magnetic 
surface stress for Cr concentrations in the range 10 % < x < 20 %. Thus, in Fe-rich 
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Figure 3. Contributions of Fe and Cr to the surface moment enhancement, Am^, 
weighted by their atomic concentration (left-hand ordinate) and magnetic surface stress 
Tmag (right-hand ordinate) . Note that Tmag equals the indicated Tmag from figure [1] 
Lines are a guide to the eye. 



Fe-Cr alloys the trend of Tmag as a function of concentration seems to be dominated by 
the magnetism of Fe. 



4-2.3. Magnetic pressure Starting from the definition of the magnetic surface stress 
(Tmag = Tnsp " Tsp) and usiug ( TTOl) . wc rcgroup all appearing terms according to 
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In the previous line, we defined the magnetic stress of the surface reference system, r™^g, 
and the magnetic stress of the bulk reference system, t^^. Notice that the definition of 
T^^ does not include any parameter related to surfaces anymore. It quantifies how much 
magnetism contributes to the bulk stress and it is thus closely related to the magnetic 
pressure known since the advent of band structure theory. [971 EH] Magnetic pressure is 
for example associated with increased atomic volumes in FM transition metals compared 
to what their volumes would be in the absence of spin-polarisation. 

The magnetic pressure of the bulk reference system per layer and per particle, 
T^g/n^, as a function of concentration for the Fe-Cr system as plotted in figure H] 



Surface parameters of ferritic iron-rich Fe-Cr alloy 

1.6 



^ 1.5 



16 



B 1.4 



M 

cs 



M 



1.3 



1.2 



1.1 



< 




^ 




' 1 ' 1 
■-■ bulk reference 




• 

• 


'^-> 


^^X 


<^-<l surface reference 

N ^\. 










^\. ^^^./\^ " 










"^^ '^\^ 










\v 'd 


1 1 






1 


, > 



0.08 







5 10 

at.% Cr X 



15 



0.07 1 

PQ 






00 



0.06 



20 



Figure 4. Correlation between the magnetic stress of the bulk reference system and 
the surface reference system per layer and per particle (left-hand ordinate) and the 
difference in the equilibrium Wigner-Seitz radius between the FM Fe-Cr system and 
the non-polarised Fe-Cr system (right-hand ordinate). Lines are a guide to the eye. 



is non- linear and tensile, indicating that without magnetism the lattice constant, or 
equivalently, the Wigner-Seitz radius would be smaller. The order of Tj^ag" ^^Y be 
connected to the difference in the Wigner-Seitz radius, r-ws? between the FM Fe-Cr 
system and the nonmagnetic Fe-Cr system. Our data shows that the equilibrium 
Wigner-Seitz radius of the nonmagnetic Fe-Cr system follows a linear concentration 
dependence, i.e, Vegard's rule (rws(Fei_a,Crj,.) ~ (1 — x)rws(F6) + a;rws(Cr)), with 
rws(Cr) > rws(Fe). As mentioned in the beginning of Sec. 14.2.11 the equilibrium 
Wigner-Seitz radius of the FM Fe-Cr system changes non-linearly as a function of the 
Cr concentration. The difference in the equilibrium Wigner-Seitz radius between the FM 
Fe-Cr alloy and the NM Fe-Cr alloy, Arws, is plotted in figure HI We find that r^^ and 
Arws are strongly correlated: t^^ as a quantitative measure for the magnetic pressure 
in the bulk Fe-Cr system correlates with the deviation of the equilibrium Wigner-Seitz 
radius from Vegard's rule. 

The magnetic pressure of the surface reference system (r™ag) tells us about the 
effect of magnetism at the surface on r compared to a magnetism-free surface. The 
surface reference systems, however, include bulk-like contributions, as effects due to 
the presence of a surface decay towards the interior of the surface reference systems 
which thus become gradually more bulk-like as the distance to the surface is increased. 
Therefore r^g and Arws are also correlated in the same way as t^^ and Arws are, see 
the plot of T^ag/""-!!! in figure HI The differences between r^^g and t^^, i.e., both the 
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absolute value and the trend as a function of x, are ascribed to both the spin-polarised 

j-surf 
mag 



and the non-spin-polarised surfaces. It is interesting to note that the maximum of r^"'"^ 



seems to be well below 5 at.% Cr while it is above 5 at.% Cr for t^^^. 

mag 

Magnetism is the driving force for the enlarged lattice parameters of the FM Fe-Cr 
system compared to the non-polarised model system. Magnetism leads on the other 
hand to increased magnetic stresses of the bulk and the surface systems. Magnetic 
stresses are compressive, i.e., have a tendency to expand the lattice. According to 
figure m T^g is larger than t^^, i.e., the magnetic contributions to r favour a larger 
surface lattice parameter than in the bulk. 

4-3. Chemically inhomogeneous Fe-Cr alloy 

There is experimental and theoretical evidence for segregation on Fe-Cr surfaces. 
According to calculations of Ropo et al for the equilibrium segregation profile of the 
(100) surface at various temperatures, the surface chemistry seems to be determined by 
the bulk configuration which leads to the peculiar threshold behaviour at approximately 
10 at.% Cr bulk content. [IHl [H] Below this value, the thermodynamically most stable 
(100) surface is essentially Fe-terminated, while it is enriched in Cr to even higher than 
bulk concentrations above this threshold up to approximately 17 at.% Cr in the bulk. 
Owning the complexity of these calculations only the surface layer concentration was 
variable. It is, however, reasonable to assume that an unconstrained computation with 
variable concentrations in a multitude of layers leads to a more complex equilibrium 
segregation profile. 

Since one nevertheless may expect that the (magnetic) contributions to both excess 
quantities originate from the topmost layers on every surface, we have good reasons 
to believe that changing only the Cr concentration of the surface layer captures the 
dominant effect on the trends of both surface energy and surface stress with surface 
alloying. It is then of course of interest to track how these surface parameters behave 
and accordingly how the stability of the surface is affected. 

The surface energies for chemically inhomogeneous concentration profiles {x^} of 
binary Ki^x^^ alloys are obtained according to [TTl [99] . 

-^surf(l^al) f^-^bulk(^) 



7({a;a}) 

(14) 



2A 



2A 

where {xa} = xi,X2, ■ ■ ■ ,Xn^ denotes the concentration of the B element within the 
layer a perpendicular to the surface, xi and x„^ being the concentrations of the two 
surface layers and Xi = Xn^-i+i due to the symmetry of the slab. The bulk effective 
chemical potential, A/ibuik(a^), equals the difference of the chemical potentials of the two 
alloy components in the bulk and is derived from the bulk energy (per atom), viz. 

dEbuik(a;) , - 

A/ibuik(a;) = — -7 • (15) 
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In order to derive the surface stress for a chemically inhomogeneous concentration 
profile, we start from ([8]) for the n-layer slab, 



-(K})= TTJ 



1 /diF/'"" - Bsi/?^." 1 



2A\ de ^ ^^, 

d (A/ibulk(x) J2lZl{Xa - X)) 



(16) 



de 



=0- 
The effective chemical potential depends on the strain. Applying a strain of the 

size e = ±0.02 led to a relative change of |A/ibuik(a;)| of approximately 10~^. Since 

Sa=i(^a — x) is at most 0.05, the second term in the previous equation containing 

the effective chemical potential is approximately by a factor of 10~^ smaller than the 

first one for the current system. Thus, we assume that the strain dependence of the 

effective chemical potential can be neglected, and r is again obtained from ( ITOj) . Thus, 

the coefficient Csurf is the one for the chemically homogeneous system with fixed bulk 

concentration x. This seems reasonable since the surface chemistry should not affect 

the bulk contributions to the surface stress. We recall that Cgurf should be identically 

zero at bulk equilibrium. 

In this work, we are mostly interested in identifying trends that arise if the 
homogeneous concentration profile is altered towards a chemically inhomogeneous 
profile. Here we were guided by the aforementioned computed equilibrium profile at 
T = K for the bcc < 100 > surfaces of Fe-Cr as reported by Ropo et al [131 [H]) i-^-, 
the concentration profiles are variable in the surface layer concentration, xi {xi = Xn^), 
and all other concentrations are fixed to the bulk value, X2 = x^ = . . . = Xn,^^^i = x. 
Two distinct bulk alloy systems, Fe95Cro5 and FegoCrio, were selected. The former has a 
predicted equilibrium surface concentration of at.% Cr, i.e., a surface Cr concentration 
lower than the bulk one, and the latter possesses a predicted equihbrium surface 
concentration of 15 at.% Cr, i.e., a surface Cr concentration larger than the bulk one. 
For the FegsCros bulk alloy system, we thus monitored 7 and r with gradual reduction 
of the Cr surface concentration from 5 to Oat.%. For the FegoCrio bulk alloy we varied 
the Cr amount at the surface between 5 and 15 at.%. 

The surface parameters for the two different bulk alloy systems are compiled in 
figure [5l For the FegsCros bulk alloy, we identify the following trend for the spin- 
polarised calculations: less surface Cr reduces both the surface energy and the surface 
stress, i.e., the Cr free surface possesses the lowest surface energy and the lowest surface 
stress. The excess surface stress (not shown) increases slightly by 0.02 J ■ m~^ when the 
surface Cr concentration is reduced from 5 at.% to Oat.%. The bulk alloys containing 
10 at.% Cr show different behaviours: the richer the surface in Cr is the smaller is the 
surface energy and the larger is the surface stress. As a result of these trends, the excess 
surface stress of the surface with xi =15 at.% Cr is by 0.03 J ■ m~^ more stable than the 
low-Cr surface. 

Our data show that an Fe-rich surface has the lowest surface energy for a bulk 
concentration below the anticipated threshold concentration, while above this threshold 
a Cr-enriched surface possesses the lowest surface energy. This finding is in qualitative 
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agreement with calculations for the surface energy of the bcc (100) facet of Fe- 
Cr from [H] done for fixed Xi = {0, 10} at.% Cr and variable bulk concentration. 
Recalling that the global trend of 7 as a function of the bulk Cr concentration shows a 
homogeneously increasing tendency in x (figure [1]) , we realise that changing the surface 
composition may alter this picture, as in the case of the FegoCrio alloy. 

We find for the investigated inhomogeneous surfaces, that a larger amount of Cr 
in the bulk and at the surface increases the surface stress r and decreases the magnetic 
surface stress Tmag (data for the inhomogeneous surface is not shown). That is, Cr 
addition drives the tendency towards smaller in-plane lattice constants at the surface 
compared to the bulk, and Cr addition has a tendency to reduce the compressive 
magnetic contribution to the total surface stress. 

5. Conclusion 

Using the EMTO method and the CPA we computed the surface energy, the surface 
stress, and the excess surface stress of the thermodynamically most stable surface facet 
(100) of the homogeneous disordered bcc Fe-Cr system in the concentration interval 
up to 20 at.% Cr. We found that the surface energy increases monotonically with Cr 
addition thereby following the rule of mixing. An increase of the bulk amount of Cr 
also increases the surface stress, which is unexpected, since the surface stress of Cr is 
considerably smaller than the one of Fe. As a result of this surprising trend, the excess 
surface stress reduces with increasing Cr concentration meaning that the (100) surface 
of Fe-Cr becomes more stable against reconstruction than the same surface of Fe. 

The reduction of the compressive magnetic contribution to the total surface stress 
(magnetic surface stress) was identified to dominate this increase of the surface stress. 
We showed further that mainly the magnetic moment enhancement of Fe is correlated 
with the behaviour of the magnetic surface stress. Thus, we conclude that mainly the 
magnetism of Fe in Fe-Cr up to 20 at.% Cr is responsible for the unexpected trend of the 
surface stress. Since the surface stress of pure Cr low-index surface facets is much smaller 
than the one of pure Fe (which was previously shown to be due to magnetism), we expect 
that Cr replaces Fe to dominate the magnetic surface stress for larger concentrations 
than 20 at.% Cr. 

We also investigated the impact of surface segregation on the surface parameters for 
the Fe95Cro5 and FegoCrio alloys. The former was previously shown to be Fe-terminated 
while the latter is enriched in Cr in vacuum. Varying only the concentration of the 
surface layer, we established the following trends: a larger amount of surface Cr increases 
the surface stress for both systems, while Cr addition raises (lowers) the surface energy 
for the bulk FegsCros (FegoCrio) alloy. 

For all investigated chemically homogeneous and inhomogeneous disordered surface 
profiles, a larger amount of Cr in the alloy favours a smaller in-plane lattice constants 
at the surface than in the bulk, and the addition of Cr shows a tendency to reduce the 
compressive magnetic contribution to the total surface stress. 



Surface parameters of ferritic iron-rich Fe-Cr alloy 21 

Acknowledgments 

The Swedish Research Council, the Swedish Steel Producers' Association, the European 
Research Council, and the Hungarian Scientific Research Fund (research project OTKA 
84078) are acknowledged for financial support. S. S. gratefully acknowledges the 
Carl Tryggers Stiftelse for Vetenskaplig Forskning and the Olle Erikssons Stiftelse for 
Materialteknik. 

References 

[1] I. Cook. Materials research for fusion energy. Nature Mat., 5:77, 2006. 

[2] R. L. Klueh and A. T. Nelson. Ferritic/martensitic steels for next-generation reactors. 

J. Nud. Mater., 371:37-52, 2007. 
[3] F. A. Garner, M. B. Toloczko, and B. H. Sencer. Comparison of swelling an irradiation 
creep behavior of fcc-austenitic and bcc-ferritic/martensitic alloys at high neutron exposure. 
J. Nucl. Mater., 276:123-142, 2000. 
[4] E. A. Little and D. A. Stow. Void-swelling in irons and ferritic steels. II. An experimental survey 

of materials irradiated in a fast reactor. J. Nucl. Mater., 87:25-39, 1979. 
[5] A. Kohyama, A. Hishinuma, D. S. Gelles, R. L. Klueh, W. Dietz, and K. Ehrlich. Low-activation 

ferritic and martensitic steels for fusion application. J. Nucl. Mater., 233-237:138-147, 1996. 
[6] G. Wranglen. An introduction to corrosion and protection of metals. Chapman and Hall, London, 

1985. 
[7] A. S. Khanna. Introduction to high temperature oxidation and corrosion. ASM International, 

Materials Park, OH, 2002. 
[8] R. N. Gunn, editor. Duplex stainless steels: Micro structure, properties and applications. Abington 

publishing, Cambridge, England, 1997. 
[9] D. E. Williams, R. C. Newman, Q. Song, and R. G. Kelly. Passivity breakdown and pitting 
corrosion of binary alloys. Nature, 350:216-219, 1991. 
[10] K. Sieradzki and R. C. Newman. A Percolation Model for Passivation in Stainless Steels. 

J. Electrochem. Soc, 133(9):1979-1980, 1986. 
[11] G. T. Burstein and P. I. Marshall. The coupled kinetics of film growth and dissolution of stainless 

steel repassivating in acid solutions. Corros. Sci, 24(5):449-462, 1984. 
[12] A. Kiejna and E. Wachowicz. Segregation of Cr impurities at bcc iron surfaces: First-principles 

calculations. Phys. Rev. B, 78:113403, 2008. 
[13] M. Ropo, K. Kokko, M. P. J. Punkkinen, S. Hogmark, J. Kollar, B. Johansson, and L. Vitos. 
Theoretical evidence of the compositional threshold behavior of FeCr surfaces. Phys. Rev. B, 
76:220401(R), 2007. 
[14] M. Ropo, K. Kokko, E. Airiskallio, M. P. J. Punkkinen, S. Hogmark, J. Kollar, B. Johansson, 
and L. Vitos. First-principles atomistic study of surfaces of Fe-rich Fe-Cr. J. Phys.: Condens. 
Matter, 23:265004, 2011. 
[15] T. P. C. Klaver, R. Drautz, and M. W. Finnis. Magnetism and thermodynamics of defect-free 

Fe-Cr alloys. Phys. Rev. B, 74:094435, 2006. 
[16] A. V. Ponomareva, E. I. Isaev, N. V. Skorodumova, Yu. Kh. Vekilov, and I. A. Abrikosov. Surface 

segregation energy in bcc Fe-rich Fe-Cr alloys. Phys. Rev. B, 75:245406, 2007. 
[17] A. V. Ruban, H. L. Skriver, and J. K. N0rskov. Surface segregation energies in transition-metal 

alloys. Phys. Rev. B, 59:15990, 1999. 
[18] M. Levesque, M. Gupta, and R. P. Gupta. Electronic origin of the anomalous segregation behavior 

of Cr in Fe-rich Fe-Cr alloys. Phys. Rev. B, 85:064111, 2012. 
[19] B. Nonas, K. Wildberger, R. Zeller, and P. H. Dederichs. Energetics of 3d Impurities on the (001) 
Surface of Iron. Phys. Rev. Lett, 80:4574, 1998. 



Surface parameters of ferritic iron-rich Fe-Cr alloy 22 

[20 

[21 



[22; 

[23; 

[24; 
[25; 

[26; 
[27; 

[28; 

[29; 

[30 
[31 

[32; 
[33; 
[34; 
[35; 

[36 

[37; 

[38; 

[39 

[4o; 

[41 

[42; 
[43; 
[44; 
[45; 



W. T. Geng. Cr segregation at the Fe-Cr surface: A first-principles GGA investigation. 

Phys. Rev. B, 68:233402, 2003. 
P. A. Dowben, M. Grunze, and D. Wright. Ssurface segregation of chromium in a Fe72Cr28(110) 

crystal Surf. Sci. Lett, 134:L524-L528, 1983. 
J. R. Lince, S. V. Didziuhs, D. K. Shuh, T. D. Durbin, and J. A. Yarmoff. Interaction of O2 

with the Feo.84Cro.i6(0Gl) surface studied by photoelectron spectroscopy. Surf. Sci., 277:43-63, 

1992. 
S. Suzuki, T. Kosaka, H. Inoue, M. Isshiki, and Y. Waseda. Effect of the surface segregation 

of chromium on oxidation of high-purity Fe-Cr aUoys at room temperature. Appl. Surf. Sci., 

103:495-502, 1996. 
C. Leygraf, G. Hultquist, S. Ekelund, and J. C. Eriksson. Surface composition studies of the (110) 

and (110) faces of monocrystalhne Feo.84Cro.i6- Surf. Sci., 46:157-176, 1974. 
G. J. Ackland. MagneticaUy Induced Immiscibihty in the Ising Model of FeCr Stainless Steel. 

Phys. Rev. Lett, 97:015502, 2006. 
P. Miillcr and A. Saiil. Elastic effects on surfaces physics. Surf. Sci. Rep., 54:157, 2004. 
H. Ibach. The role of surface stress in reconstruction, epitaxial growth and stabilization of 

mesoscopic structures. Surf. Sci. Rep., 29:193, 1997. 
H. Ibaeh. Erratum: The role of surface stress in reconstruction, epitaxial growth and stabilization 

of mesoscopic structures. Surf. Sci. Rep., 35:71, 1999. 
W. Haiss. Surface stress of clean and adsorbate-covered solids. Rep. Prog. Phys., 64(5):591, 2001. 
G. Wulff. Zur Frage der Geschwindigkeit des Wachstums und der Auflosung der Krystallflachen. 

Z. Kristall. Mineral, 34:449-530, 1901. 
C. E. Bach, M. Giesen, H. Ibaeh, and T. L. Einstein. Stress relief in reconstruction. 

Phys. Rev. Lett, 78:4225-4228, 1997. 
S. Olivier, A. Saul, and G. Treglia. Relation between surface stress and (1 x 2) reconstruction for 

(110) fee transition metal surfaces. Appl. Surf. Sci., 212-213:866-871, 2003. 
P. Wynblatt and R. C. Ku. Surface energy and solute energy effects in surface segregation. 

Surf. Sci., 65:511-531, 1977. 
H. Ibach. The relation between the strain-dependence of the heat of adsorption and the coverage 

dependence of the adsorbate induced surface stress. Surf. Sci., 556:71-77, 2004. 
M. Alden, H. L. Skriver, S. Mirbt, and B. Johansson. Calculated surface-energy anomaly in the 

3d metals. Phys. Rev. Lett, 69:2296, 1992. 
M. P. J. Punkkinen, S. K. Kwon, J. KoUar, B. Johansson, and L. Vitos. Compressive surface 

stress in magnetic transition metals. Phys. Rev. Lett., 106:057202, 2011. 
M. P. J. Punkkinen, Q. M. Hu, S. K. Kwon, B. Johansson, J. Kollar, and L. Vitos. Surface 

properties of 3d transition metals. Philosophical Magazine, 91:3627, 2011. 
M. Blanco- Rey and S. J. Jenkins. Surface stress in d-band metal surfaces. J. Phys.: Condens. Mat- 
ter, 22:135007, 2010. 
T. Ossowski and A. Kiejna. Density functional study of surface properties of chromium. Surf. 

Set., 602:517, 2008. 
R. Lamber, S. Wetjen, and N. I. Jaeger. Size dependence of the lattice parameter of small palladium 

particles. Phys. Rev. B, 51:10968, 1995. 
H. J. Wasserman and J. S. Vermaak. On the determination of a lattice contraction in very small 

silver particles. Surf Set., 22:164-172, 1970. 
H. J. Wasserman and J. S. Vermaak. On the determination of the surface stress of copper and 

platinum. Surf Sci, 32:168-174, 1972. 
C. SoUiard and M. Flueli. Surface stress and size effect on the lattice parameter in small particles 

of gold and platinum. Surf Sci., 156:487-494, 1985. 
C. W. Mays, J. S. Vermaak, and D. Kuhlmann-Wilsdorf. On surface stress and surface tension. 

II. Determination of the surface stress of gold. Surf. Sci., 12:134-140, 1968. 
J. J. Metois and P. Miillcr. Absolute surface energy determination. Surf. Sci., 548:13-21, 2004. 



Surface parameters of ferritic iron-rich Fe-Cr alloy 23 



[46; 
[47; 

[48; 

[49; 

[50 

[51 
[52; 
[53; 
[54; 
[55; 
[56 
[57; 

[58; 



[59 
[60 

[61 

[62; 

[63 

[64; 

[65 



[66 

[67 
[68; 
[69 

[7o; 

[71 

[72; 



R. Shuttleworth. The Surface tension of solids. Proc. Phys. Soc. A, 63:444, 1950. 

R. J. Needs. Calculations of the surface stress tensor at aluminum (111) and (110) surfaces. 

Phys. Rev. Lett, 58:53, 1987. 
J. C. Boettger. Nonconvergence of surface energies obtained from thin-film calculations. 

Phys. Rev. B, 49:16798, 1994. 
V. Fiorcntini and M. Mcthfessel. Extracting convergent surface energies from slab calculations. 

J. Phys.: Condens. Matter, 8:6525, 1996. 
J. Kollar, L. Vitos, J. M. Osorio-Guillen, and R. Ahuja. Calculation of surface stress for fee 

transition metals. Phys. Rev. B, 68:245417, 2003. 
K. Kadas, Z. Nabi, S. K. Kwon, L. Vitos, R. Ahuja, B. Johansson, and J. Kollar. Surface relaxation 

and surface stress of Ad transition metals. Surf. Sci., 600:395, 2006. 
S. K. Kwon, Z. Nabi, K. Kadas, L. Vitos, J. Kollar, B. Johansson, and R. Ahuja. Surface energy 

and stress release by layer relaxation. Phys. Rev. B, 72:235423, 2005. 
R. C. Cammarata. Continuum model for surface reconstruction in (111) annd (100) oriented 

surfaces of fee metals. Surf. Sci., 279:341-348, 1992. 
V. Zolyomi, L. Vitos, S. K. Kwon, and J. Kollar. Surface relaxation and stress for 5d transition 

metals. J. Phys.: Condens. Matter, 21:095007, 2009. 
L. Vitos. Computational quantum mechanics for materials engineers: the EMTO method and 

applications. Springer- Verlag, London, 2007. 
L. Vitos. Total-energy method based on the exact muffin-tin orbitals method. Phys. Rev. B, 

64:014107, 2001. 
L. Vitos, H. L. Skriver, B. Johansson, and J. Kollar. Application of the exact muffin-tin orbitals 

theory: the sperical cell approximation. Comp. Mat. Sci, 18:24, 2000. 
J. Zabloudil, R. Hammerling, L. Szunyogh, and P. Weinberger. Electron scattering in solid matter: 

A theoretical and computational treatise, volume 147 of Springer series in solid-state sciences. 

Springer- Verlag, 2005. 
J. Korringa. On the calculation of the energy of a Bloch wave in a metal. Physica, 13:392, 1947. 
W. Kohn and N. Rostoker. Solution of the Schrodinger equation in periodic lattices with an 

application to metallic lithium. Phys. Rev., 94:1111, 1954. 
P. Soven. Coherent-potential model of substitutional disordered alloys. Phys. Rev., 156:809, 1967. 
B. L. Gyorffy. Coherent-potential approximation for a nonoverlapping-mufhn-tin-potential model 

of random substitutional alloys. Phys. Rev. B, 5:2382, 1972. 
J. S. Faulkner. The modern theory of alloys. Prog. Mater. Sci., 27:1-187, 1982. 
M. Zwierzycki and O. K. Andersen. The overlapping muffin-tin approximation. Act. Phys. Pol. A, 

115:64, 2009. 
J. Kollar, L. Vitos, and H. L. Skriver. Electronic structure and physical properties of solids: the 

uses of the LMTO method, volume 535 of Lecture notes in physics. Springer- Verlag, Berlin, 

2000. 
L. Vitos, J. Kollar, and H. L. Skriver. Full charge-density calculation of the surface energy of 

metals. Phys. Rev. B, 49:16694, 1994. 
L. Vitos, J. Kollar, and H. L. Skriver. Ab initio full charge-density study of the atomic volume of 

a-phase Fr, Ra, Ac, Th, Pa, U, Np, and Pu. Phys. Rev. B, 55:4947, 1997. 
L. Vitos, J. Kollar, and H. L. Skriver. Full charge-density scheme with a kinetic-energy correction: 

Application to ground-state properties of the Ad metals. Phys. Rev. B, 55:13521, 1997. 
L. Vitos, I. A. Abrikosov, and B. Johansson. Anisotropic lattice distortions in random alloys from 

first-principles theory. Phys. Rev. Lett., 87:156401, 2001. 
A. V. Ruban, P. A. Korzhavyi, and B. Johansson. First-principles theory of magnetically driven 

anomalous ordering in bee Fe-Cr alloys. Phys. Rev. B, 77:094436, 2008. 
V. I. Razumovskiy, A. V. Ruban, and P. A. Korzhavyi. First-principles study of elastic properties 

of Cr- and Fe-rich Fe-Cr alloys. Phys. Rev. B, 84:024106, 2011. 
P. Olsson, I. A. Abrikosov, L. Vitos, and J. Wallenius. Ab initio formation energies of Fe-Cr alloys. 



Surface parameters of ferritic iron-rich Fe-Cr alloy 24 



[73; 
[74; 

[75; 
[76; 

[77 

[78; 

[79; 
[80; 

[81 

[82; 
[83; 
[84; 
[85; 
[86; 

[87 



[90 
[91 

[92; 

[93 

[94 
[95 



J. Nucl. Mater., 321:84, 2003. 
P. Olsson, I. A. Abrikosov, and J. Wallenius. Electronic origin of the anomalous stability of Fe-rich 

bcc Fe-Cr alloys. Phys. Rev. B, 73:104416, 2006. 
H. Zhang, B. Johansson, and L. Vitos. An ab initio calculations of elastic properties of bcc Fe-Mg 

and Fe-Cr random alloys. Phys. Rev. B, 79:224201, 2009. 
P. A. Korzhavyi, A. V. Ruban, J. Odqvist, J.-O. Nilsson, and B. Johansson. Electronic structure 

and effective chemical and magnetic exchange interactions in bcc Fe-Cr alloys. Phys. Rev. B, 

79:054202, 2009. 
P. A. Korzhavyi, A. V. Ruban, I. A. Abrikosov, and H. L. Skriver. Madelung energy for random 

metallic alloys in the coherent potential approximation. Phys. Rev. B, 51:5773, 1995. 
A. V. Ruban and H. L. Skriver. Screened Coulomb interactions in metallic alloys. I. Universal 

screening in the atomic-sphere approximation. Phys. Rev. B, 66:024201, 2002. 
A. V. Ruban, S. I. Simak, P. A. Korzhavyi, and H. L. Skriver. Screened Coulomb interactions 

in metallic alloys. II. Screening beyond the single-site and atomic-sphere approximation. 

Phys. Rev. B, 66:024202, 2002. 
I. Mirebeau, M. Hennion, and G. Parette. First measurement of short-range-order inversion as a 

function of concentration in a transition alloy. Phys. Rev. Lett., 53:687, 1984. 
I. Mirebeau and G. Parette. Neutron study of the short range order inversion in Fci-ajCr^:. 

Phys. Rev. B, 82:104203, 2010. 
J. P. Perdew, K. Burke, and M. Ernzerhof. Generalized gradient approximation made simple. 

Phys. Rev. Lett, 77:3865, 1996. 
J. P. Perdew, K. Burke, and M. Ernzerhof. Erratum: Generalized gradient approximation made 

simple. Phys. Rev. Lett., 78:1396, 1997. 
J. P. Perdew and Y. Wang. Accurate and simple analytic representation of the electron-gas 

correlation energy. Phys. Rev. B, 45:13244, 1992. 
M. Asato, A. Settcls, T. Hoshino, T. Asada, S. Bliigel, R. Zcller, and P. H. Dederichs. Full-potential 

KKR calculations for metals and semiconductors. Phys. Rev. B, 60:5202, 1999. 
K. Kocpcrnik and H. Eschrig. Full-potential nonorthogonal local-orbital minimum-basis band- 
structure scheme. Phys. Rev. B, 59:1743, 1999. 
P. R. Watson, M. A. van Hove, and K. Hermann. Atlas of surface structures, volume lA and IB. 

Am. Chem. Soc./AIP, 1994. ; NIST Standard Ref. Database 42 (1995). 
P. Blonski and A. Kiejna. Structural, electronic, and magnetic properties of bcc iron surfaces. 

Surf Sci., 601:123-133, 2007. 
G. Chiarotti, editor. The Landolt-Bornstein Database, chapter 2.2.2.1 Metals. SpringerMaterials, 

2012. Tables 7-10, Figs. 1-7. 
Z. Q. Wang, Y. S. Li, F. Jona, and P. M. Marcus. Epitaxial growth of body-centered-cubic nickel 

on iron. Solid State Commun., 61(10):623-626, 1987. 
P. Villars and L. D. Calvert. Pearson's handbook of crystallographic data for intermetallic phases, 

volume 4. ASM international. Materials Park, Ohio, 2nd edition, 1991. 
H. P. J. Wijn, editor. The Landolt-Bornstein Database, volume 32A of The Landolt-Bornstein 

Database, chapter 1.1.2.7 Magnetic moments, Compton profiles. SpringerMaterials, 1997. 
J. P. Perdew, J. A. Chcvary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, and 

C. Fiolhais. Atoms, molecules, solids, and surfaces: applications of the generalized gradient 

approximation for exchange and correlation. Phys. Rev. B, 46:6671, 1992. 
J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, and 

C. Fiolhais. Erratum: Atoms, molecules, solids, and surfaces: applications of the generalized 

gradient approximation for exchange and correlation. Phys. Rev. B, 48:4978, 1993. 
C. J. Fall, N. Binggeli, and A. Baldereschi. Anomaly in the anisotropy of the aluminum work 

function. Phys. Rev. B, 58:R7544, 1998. 
W. R. Tyson and W. A. Miller. Surface free energies of solid metals: estimation from liquid surface 

tension measurements. Surf. Sci., 62:267, 1977. 



Surface parameters of ferritic iron-rich Fe-Cr alloy 25 

[96] M. J. S. Spencer, A. Hung, I. K. Snook, and I. Yarovsky. Density functional theory study of the 

relaxation ans energy of iron surfaces. Surf. Sci., 513:389-398, 2002. 
[97] O. K. Andersen, O. Jepsen, and D. Glotzel. Highlights of condensed-matter physics, chapter 

Canonical description of the band structure of metals, pages 59-176. Proceedings of the 

international school of physics 'Enrico Fermi'. North-Holland, Bassani, F., Fumi, F., and Tossi, 

M. P., Eds., Amsterdam, 1985. 
[98] O. K. Andersen, J. Madsen, U. K. Poulsen, O. Jepsen, and J. KoUar. Magnetic ground state 

properties of transition metals. Physica, 86-88B:249-256, 1977. 
[99] L. V. Pourovskii, A. V. Ruban, B. Johansson, and I. A. Abrikosov. Antisite-defect-induced surface 

segregation in ordered NiPt alloy. Phys. Rev. Lett, 90:026105, 2003. 



